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Using Monte Carlo techniques, Ising models with ferro- 
magnetic nearest-neighbor interactions on a simple cubic lat- 
tice are studied. At the surface transition, the critical ex- 
ponent /?2 of the edge magnetization is found to be non- 
universal, depending on the edge and edge-surface couplings, 
in contrast to the situation at the ordinary transition. Results 
are compared to those for two-dimensional Ising magnets with 
chain and ladder defects. 

Pacs numbers: 05.50+q, 68.35.Rh, 75.40.Mg 



I. Introduction 

Critical phenomena occur not only in the bulk of a sys- 
tem, but also at its surfaces and edges. To be specific, 
let us consider the magnetization as the order parameter. 
Then, there are two typical scenarios: (a) bulk, mb, sur- 
face, mi, and edge, 171%, magnetizations may order at the 
same temperature ('ordinary transition'), but with dif- 
ferent power-laws, and (b) surface and edge may order 
first ('surface transition') due to strong surface couplings, 
followed by ordering of the bulk magnetization ('extraor- 
dinary transition') at a lower temperature. Surface sin- 
gularities at ordinary and surface transitions have bee 
studied extensively for various systems, theoretically^ 
as well as experimentally.^ Edge critical behavior, on the 
other hand, has attracted less attention. 

Some years ago, Cardy noted and calculated the de- 
pendence of edge critical exponents on the opening angle 
between the surfaces forming the edge, using mean field 
theory and renormalization group theory of fiijst order 
in e for 0(n) models at the ordinary transitions Subse- 
quent work dealt mainly with edge criticality at the or- 
dinary transition as well, including high temperature se- 
ries expansionsEl, exact analyses of two-dimensional Ising 
modelsa, and applications to polymers.Q 

In recent Monte Carlo simulations of three-dimensional 
Ising models, we confirmed and refined the estimates ob- 
tained from renormalization group calculations and high 
temperature series expansions on edge critical exponents 
at the ordinary transition, especially on 02 describing the 
vanishing .of the edge magnetization mi on approach to 
criticality.a The value of /?2 varies with the opening an- 
gle, but not with the bulk, surface and edge coupling 
constants which are assumed to be ferromagnetic and of 
short range. 

Edge critical properties at the surface transition have 
been largely overlooked. Actually, we are aware of merely 
one exception, where they have been mentioned briefly 
and qualitatively.EI However, that case deserves to be in- 



vestigated, in particular, by Monte Carlo techniques in 
the framework of Ising models, too. At the surface tran- 
sition in three-dimensional systems, the critical fluctu- 
ations are essentially two-dimensional, and the surface 
critical exponents reflect the reduced dimensionalityJj'EI 
On the other hand, the edge presents a local pertur- 
bation, acting presumably like the much studied chain 
or ladder defect sEJfli-1 in two-dimensional Ising systems. 
Following those studies, intriguing non-universal edge 
critical behavior, depending on the interactions at the 
edge, may be expected, in contrast to the situation at the 
ordinary transition. Furthermore, previous Monte Carlo 
simulations on the critical jjuxinerties of Ising models at 
surfaces, edges, and corneraali3ii3 show that present sim- 
ulational techniques, especially cluster-flip algorithms, 
allow to estimate reliably edge critical properties. Ac- 
cordingly, the simulational results may provide guidance 
for future analytical investigations. Finally, the find- 
ings may encourage experimental search for suitable mag- 
nets or alloysu, where the surface coupling is sufficiently 
strongly enhanced compared to the bulk coupling. 

The outline of the article is as follows. In the next 
section, the model and its phase diagram as well as the 
Monte Carlo method will be introduced. In section 3, 
we shall discuss magnetization profiles and the critical 
exponent of the edge magnetization at the surface tran- 
sition in the three-dimensional Ising model for various 
couplings near the edges. A brief summary concludes 
the paper. 

II. Ising model with surfaces and edges 

We study nearest-neighbor Ising models on simple cu- 
bic lattices with ferromagnetic interactions. The Hamil- 
tonian may be written in the form 

'ri. ^ JbSxyzSx'y' z' ^ J s^xyz^x 1 y 1 z 1 

bulk surface 

^ Jes^xyzSx'y' z' ^ J e^xyz^x' y' z' (1) 

edge— surf ace edge 

where the sums run over bonds between neighboring 
spins, Sxyz — ±1) with coupling constants to be specified 
below. Free boundary conditions hold for the spins in the 
yz- and xz-surface planes, while the spins in the first and 
last xy-planes are connected by periodic boundary con- 
ditions, see Fig.l. Thereby the edges formed by the inter- 
secting free surface planes are oriented along the z-axis. 
The total number of sites or spins is L x M x N, where 
L, M, and N correspond to the x, y,and z-directions, re- 
spectively. In the following, we assume L — M, i.e. there 
are L 2 spins in each plane perpendicular to the edges. 
The pairs of spins in the Hamiltonian (1) are situated 
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either on edges with the edge coupling J e , on edge and 
surface sites coupled by J es , on surface sites with the 
interaction J s , or with at least one of the spins in the 
interior of the system interacting with the bulk coupling 
Jb, as depicted in Fig. 1. 

The magnetization per site for lines parallel to the z- 
axis, mi(x,y), is defined as 

m l {x,y)=<\Y d S xvz \> /N (2) 

z 

summing over all spins of a line, with x and y being fixed. 
The brackets denote thermal averages. The absolute 
value is taken to avoid vanishing magnetizations for finite 
systems, as usual. The edge magnetization, ma, is then 
identical to m;(l,l) = mi(l,L) — mi(L,l) = mi(L,L). 
The surface magnetization, mi, is given by m;(l,Lc) = 
mi(L c ,L) = mi(L c ,l) = ™>l(L,Lc), where L c refers 
to the one or two center lines, with Lc = (L + l)/2 or 
Lc = L/2 ± 1. The profile of the line magnetization at 
the surface, to/ s (i), is described by 

mu(i) = (mi(x,L) + mi(x, 1) + m ; (l,y) + m;(L,y))/4 

(3) 

with, setting the lattice constant equal to one, i — x — 
y = 1, 2, L (certainly, the profile is symmetric around 
i = Lc)- The bulk magnetization, nib, corresponds to 
mi(Lc,L c )- 

Of course, the line magnetizations may be affected 
by finite size effects. The thermodynamic limit is ap- 
proached when L, M, N — > oo. 

To elucidate bulk, surface, and edge properties of the 
model, we considered a few additional quantities. In par- 
ticular, we calculated the energy of the surface spins, 
E s , the energy of the edge spins, E e , the edge sus- 
ceptibility X2*i (defined as the fluctuation of the edge 
magnetization!]), and higher moments of various mag- 
netizations (Binder cumulants). 

The linear dimensions, L and AT, of the L 2 x N sys- 
tems ranged from L — 10 to L — 80, and from N = 10 to 
N = 160. We used the efficient one-cluster-flip Monte 
Carlo algorithm, which reduces significantly critical slow- 
ing down. Typically, averages were taken over a few 10 4 
clusters, after equilibration. Error bars resulted from 
sampling over several realizations. 

The phase diagram of the Ising model (1) is deter- 
mined by tigpuatio of the surface to the bulk coupling, 
r = J s /J b mHt5 see Fig. 2. At r < r c » 1.50, the 
model displays an ordinary transition, with bulk, to;,, and 
surface magnetizations, mi, ordering simultaneously at 
the bulk transition temperature, T c = 4.5115... J/fcsO : E!l 
The surface critical exponents are universal, i.e. inde- 
pendent of the ratio r. For instance, the surface mag- 
netization toi vanishes on approach to T c as toi oc t 131 , 
where t = \T — T C U^T C is the reduced temperature, with 
0i = 0.80 ± 0.0lli3 (this value is robust against ran- 
domness in the surface couplings and corrugations of the 



surfaceuIrEil) . At r = r c « 1.50, bulk and surface still 
order at the same temperature, T c , but the surface crit- 
ical exponents at this 'special PftWt' are in a different 
universality class, e.g., Pi ss 0.24.ErH 

At J s > r c Jb, the surface magnetization mi orders 
at the surface transition, with the critical temperature 
T s being higher than that of the bulk ordering which 
occurs at the extraordinary transition, T c . At T s , the 
surface critical fluctuations are of two-dimensional char- 
acter, and the surface critical exponents are those of 
the two-dir&ensional Ising model with nearest-neighbor 
interactionsoa, i.e., for example, @\ — 1/8. 

Because edges are one-dimensional, and all couplings 
in the models are of short range, one expects edge quan- 
tities to become singular at the surface transition, T a , 
but the edge critical exponents may differ from the cor- 
responding surface valuesa (similar to the critical behav- 
ior at onp-dimensional surfaces in two-dimensional Ising 
modelstj) . Most of the Monte Carlo simulations near T s 
were done at fixed ratio r = J s /Jb = 2, varying then the 
nearest neighbor interactions along the edges, J e , and be- 
tween edge and surface spins, J es , with J e ranging from 
to 2J S , and J es ranging from 0.5J S to 2J S . 

To determine edge critical exponents accurately, T s 
needs to be dctcraiincd accurately. Using standard 
finite-size analyses^ (with the scaling variable Nt, im- 
plying that the surface fluctuations are governed by the 
surface correlation length which diverges like T s 
may be obtained from the location of turning points or 
extrema in various thermal quantities. As illustrated and 
indicated in Fig. 3, different quantities give, indeed, con- 
sistent estimates. For r = 2, we find T s = 4.9575±0.0075. 
Note that other critical temperatures of the surface tran- 
sition, needed to map the phase diagram, Fig. 2, have 
been obtained with less accuracy, by doing finite-size 
analyses on the location of the turning point of the sur- 
face magnetization for systems of moderate sizes, with 
up to 40 3 sites. 

III. Magnetization profiles and edge critical ex- 
ponents 

The profiles of the line magnetization at the surface, 
mi s (i), equation (3), reflect the influence of bulk spins 
close to the surfaces and edges as well as the strength of 
the couplings near the edges, J e and J es . 

Typical profiles are shown in Fig. 4. At J e = J es = 
J s {= 2Jb), i.e., for equal edge, edge-surface, and surface 
couplings, the effect of the bulk spins near the surfaces 
and edges on the temperature dependence of the profile 
is illustrated in Fig. 4a. At low temperatures, mi s (i) 
increases monotonically with i, i.e. the distance from the 
edge. The magnetization of the surface spins, which are 
connected directly to an ordered bulk spin, is enhanced 
compared to the edge magnetization, TO2 = m; s (l). Edge 
spins have the coordination number four, while it is five 
for the surface spins. In contrast, on approach to T s , 
the ordering of the spins falls off quickly by going from 
the surface to the bulk, and the surface magnetization, 
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mi = mi s (Lc), is pulled down below the edge magne- 
tization by the coupling to the bulk spin. Roughly at 
temperatures above T c , an interesting non-monotonic be- 
havior in the profile shows up, with a maximum close 
to the edge. It may be explained by the fact that sur- 
face spins next to, but on different sides of an edge, are 
more strongly connected to each other, through the same 
neighboring bulk spin, than spins on the flat part of the 
surface with the separation distance of two. 

The profile approaches the surface magnetization mi 
nearly exponentially already at rather small distances 
from the maximum, i.e. \mi s (i) — m\\ oc exp(— ai)\ a 
is expected to be, as T — > T s , the inverse correlation 
length.™ 

In Fig. 4b, the influence of the edge and edge-surface 
couplings, J e and J es , on the profile, at fixed tempera- 
ture, is illustrated. The trends may be readily under- 
stood. By increasing the coupling, J e , along an edge, the 
edge magnetization will rise. Likewise, a larger interac- 
tion, J es , of the edge spins with the neighboring surface 
spins will lead to an increase of the edge magnetization. 
In that case, the non-monotonic shape of the profile close 
to the edge sets in at lower temperatures. 

Near the surface transition, where the critical fluctua- 
tions are of two-dimensional character, the edge acts like 
a defect line in an essentially two-dimensional bulk Ising 
model; J e corresponds to a 'chain defect' and J es to a 
'ladder-type defect'.Ej The change in the topology at the 
edge compared to the surface amounts to a complicated 
ladder- type defect. 

This analogy will be crucial in explaining the simula- 
tional findings on the critical exponent (32 of the edge 
magnetization, defined by 
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(4) 



where t is the reduced temperature of the surface tran- 
sition, t = \T S — T\/T s , t — xJk To estimate $2, we 
consider the effective exponentlaO 



Peff = dhiT^/dlnf 



(5) 



In analysing Monte Carlo data, the derivative is replaced 
by a difference at discrete temperatures. (3 e ff approaches 
/?2 as t — > 0, provided finite-size effects can be ne- 
glected. 

To check whether the edge magnetization 7712 is af- 
fected by the finite size of the Monte Carlo system, one 
may proceed as follows. First, one chooses, at fixed tem- 
perature, the linear dimensions L and N to be sufficiently 
large to reproduce the numerically to a high degree of ac- 
curacy known thermodynamic bulk, mi, = mi{Lc,Lc), 
and surface, m\ — mi s (Lc), magnetizations. In a sec- 
ond step, TV may be enlarged to search for further size 
dependences in m,2a 

In Fig. 5, the effective exponent (3 e ff of the edge mag- 
netization is displayed for various system sizes, in the case 
of equal surface and edge couplings J e = J es = J s = Ub- 



Obviously, (3 e ff increases on lowering the reduced tem- 
perature before it acquires, at t < 0.15, nearly a plateau. 
Finally, on further approach to T s , it decreases to zero 
due to a non-vanishing edge magnetization at the criti- 
cal point in finite systems. Monte Carlo data which are 
largely free of finite-size dependences have been obtained 
down to t w 0.06, using systems with up to 80 3 spins. 
Assuming the plateau-like behavior, 0.06 < t < 0.15, 
to hold even closer to T s (which would imply that cor- 
rections to scaling, i.e. to the asymptotic power-law, 
equation (4), are small), we arrive at the estimate of the 
critical exponent $2 for this set of coupling constants, 
f3 2 = 0.095 ± 0.005. Error bars in Fig. 5 include both 
sample averaging and the uncertainty in T s , see above. 
Here and in the following, the error bar for the value of /?2 
is rather subjective, based on 'reasonable' extrapolation 
of the data, which are close to those in the thermody- 
namic limit, to t — > 0. 

The critical exponent P2 is significantly lower than the 
critical exponent of the surface magnetization, (3i = 1/8, 
reflecting the fact that the edge magnetization is, on ap- 
proach to T s , larger than the surface magnetization. We 
may relate this finding to results for two-dimensional 
nearest-neighbor Ising model with a ladder defect. In 
that model, the couplings are identical, say, J > 0, 
throughout the system with the exception of the cou- 
plings, Ji, of the spins in one row to the spins in a neigh- 
boring row. The ladder rows correspond to the edge, and 
the ladder couplings, J/, reflect not only the edge-surface 
interactions J es , but also the change in topology (or con- 
nectedness to bulk spins) at the edge compared to the 
surface. For the two-dimensional Ising model with a lad- 
der defect, the critical exponent (3i of the magnetization 
in the ladder rows has been determined exactlytlrEj, 



(3i = 2arctan 2 (K ; _1 )/7r 2 



(6) 



with m = tanh(J ; /(A: B T2 ( i))/tanh(J/(fc B T2 ( i)), where 
T2d is the transition temperature. Thence, the critical ex- 
ponent of the ladder row magnetization varies monoton- 
ically with the strength of the coupling J\ , being smaller 
than 1 /8 when J/ > J. Comparing the critical exponents 
[3i and /?2, one may attribute an effective ladder cou- 
pling Ji to the edge. Our case, J e = J es = J s = (2J&), 
corresponds to Ji > J(— J s ). Of course, the non- 
monotonic profiles of the line magnetization suggest that 
a closer analogy between edge properties and descriptions 
by two-dimensional Ising models with defect lines would 
require more complicated, extended ladder-type-defects, 
which have not been treated analytically so far.tfl 

At any rate, exact results on two-dim p.yi sinnal Ising 
models with defect lines of different typesoliilllj provide 
a useful framework to discuss our simulational findings on 
the dependence of P2 on the edge and edge-surface cou- 
plings, J e and J es . Obviously, J es is intimately related 
to Ji , while changing J e may be interpreted as modifying 
nearest-neighbor couplings, J c h, along a single defect line 
or chain in the two-dimensional Ising model. For a chain 
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defect, the critical exponent, /3 c h, of theujaiagnetization 
in the defect chain varies continuously asEilO 

f3 ch = 2arctan 2 (K c/l )/vr 2 (7) 

with K ch = exp[-2 (J ch / \k B T 2d ) - J/{k B T 2d ))\. Again, 
the critical exponent changes monotonically, being larger 
than 1/8 when J c h < J. 

Pertinent results on the edge critical exponent (3 2 are 
shown in Figs. 6 and 7. In Fig. 6, the effect of weakening 
the interactions between the edge and the surface, J es , is 
demonstrated. As expected from equation (6), lowering 
that coupling leads to an increase in the value of j3 2 . 
Only Monte Carlo data are depicted which seem to be 
unaffected by finite-size dependences, simulating systems 
with up to 80 3 sites. As before, the effective exponent 
(3 e ff exhibits a broad plateau in t, which may allow one 
to extrapolate the asymptotic exponent quite accurately. 
Our estimates, for J e = J s = 2Jt, are (3 2 = 0.176 ± 0.005 
at J es = J e /2 and /3 2 = 0.244 ± 0.005 at J es = J e /4. 

Similarly, weakening of the edge interaction J e tends to 
increase the critical exponent of the edge magnetization, 
see equation (7). In Fig. 7, the effective exponent /3 e // is 
shown, using again only those simulational data for the 
edge magnetization which approximate closely the infi- 
nite system. In this case, the effective exponent depends 
on approach to T s nearly linearly on the reduced tem- 
perature t, signalling rather strong corrections to scaling 
for m 2 . By extrapolating the findings for (3 e ff to the 
critical point, t — > 0, we estimate, for J es — J s = 2Jb, 
/3 2 = 0.170 ± 0.005 in the limiting case J e = 0, and 
02 = 0.127 ± 0.005 for J e = J es /2. The close agreement 
with the value of the surface critical exponent /3i = 1/8, 
in the latter case, is rather fortuitous, due to a compen- 
sation of a reduction in /3 2 following from the increase 
in the effective ladder coupling stemming from the edge 
topology (as discussed above), and of an enhancement 
in p 2 following from the decrease in the strength of the 
chain coupling. 

The non-universal character of the edge transition 
driven by the surface transition is also confirmed by our 
findings for the critical properties of the edge susceptibil- 
ity X22- However, because the estimates are quite rough, 
we refrain from quoting specific values. 

The Binder cumulant at T s for the edge magnetiza- 
tion depends on the edge and edge-surface couplings as 
well. However, this may not be interpreted as addi- 
tional evidence for non-universality, because that cumu- 
lant changes with, e.g., the ratio of the surface to bulk 
interactions r = J s /Jb even at the ordinary tcarisition, 
while the critical exponents are universal there.E-3 

IV. Summary 

We studied magnetization profiles and the critical ex- 
ponent (3 2 of the edge magnetization of three-dimensional 
nearest-neighbor Ising models near the surface transi- 
tion, T s (with the surface coupling J s being twice as 
large as the bulk coupling). We did large-scale simu- 



lations, using the single-cluster-flip Monte Carlo algo- 
rithm. In particular, the influence of the edge, J e , and 
edge-surface, J es , couplings has been investigated. 

The profiles of the line magnetization at the sur- 
face suggest that the edge acts like a rather compli- 
cated, extended ladder-type defect in an essentially two- 
dimensional Ising magnet. The reduced dimensionality 
is due to the two-dimensional character of the critical 
fluctuations at the surface transition. 

The edge critical exponent f3 2 varies continuously with 
J e and J es . If J e and J es are equal to the surface cou- 
pling, (3 2 (=0.095 ± 0.005) is significantly lower than the 
critical exponent of the surface magnetization, (3\ = 1/8. 
When weakening the couplings at the edge, f3 2 increases 
monotonically. Eventually, [3 2 will be clearly larger than 
1/8, in accordance with related exact results for two- 
dimensional Ising models with defect lines. The decrease 
of J es has a more pronounced impact on (3 2 ; the value 
of the edge critical exponent may exceed appreciably the 
limiting value obtained at vanishing J e . By weakening 
J e , strong corrections to scaling for the edge magnetiza- 
tion are observed. 

It may be of interest to check the present Monte Carlo 
results by applying other methods, such as renormaliza- 
tion group techniques. Likewise, it may be feasible to 
mimic the edge by a realistic ladder-type defect in a, 
possibly, exactly solvable two-dimensional Ising model. 
Finally, the non-universality of the critical behavior is 
expected to show up in additional dependences of the 
critical exponents on details of, for instance, the lattice 
structure and the range of interactions. 
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FIG. 1. Sketch ol the lattice structure and the coupling 
constants at the edges, J e and J es , at the surface, J s , and in 
the bulk, Jt- 
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FIG. 2. Phase diagram of the three-dimensional Ising 
model, in the (fcsT/Jt, J s / Jb)-plane. The dashed curve de- 
notes the surface transition line, T s . 
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FIG. 3. Monte Carlo data on temperature dependence of 
the edge susceptibility, X22, and of the temperature deriva- 
tives of the edge energy, E e , surface energy, E 3 , surface mag- 
netization, mi, as well as edge magnetization, 7712, near T s , 
at J e — Jes = Js = 2Jb for a system of size 40 2 x 160. 
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FIG. 4. Profiles of the line magnetization mi s (i) at (a) 
Je — Jes — J s — 2,h, with various temperatures, and (b) 
at fixed temperature, knT/Jt, — 4.8, with various edge and 
edge-surface couplings (J s — 2Jb). Ising models of size 40 3 
have been simulated. 



FIG. 5. Effective exponent /3 e // of the edge magnetiza- 
tion as a function of reduced temperature t = \T S — T\/T s 
at J e — J es = J s = 2J b for Monte Carlo systems of various 
sizes. 
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FIG. 6. Effective exponent /3 e // versus reduced tem- 
perature t for different edge-surface couplings J es , with 
Je = Js = 2Jt. System sizes have been adjusted to circum- 
vent finite-size dependences, with up to 80 3 spins. 
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FIG. 7. Effective exponent /3 e // versus reduced tempera- 
ture t for different edge couplings J e , with J es — J s = 2Jb- 
System sizes for the simulations have been chosen to avoid 
finite-size effects, with up to 80 3 spins close to the surface 
transition. 
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